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Time-dependent Hartree-Fock theory is used to describe 
density oscillations of symmetry-unrestricted two-dimensional 
nanostructures. In the small amplitude limit the results repro- 
duce those obtained within a perturbative approach such as 
the linearized time-dependent Hartree-Fock one. The nonlin- 
ear regime is explored by studying large amplitude oscillations 
in a non-parabolic potential, which are shown to introduce a 
strong coupling with internal degrees of freedom. This exci- 
tation of internal modes, mainly of monopole and quadrupole 
character, results in sizeable modifications of the dipole ab- 
sorption. 

PACS 73.20.Dx, 72.15.Rn 



Hartree-Fock (HF) theory is one of the most success- 
ful microscopic schemes to address finite Fermi systems, 
and usually the starting point of more refined approx- 
imations. Shortly after its proposal Dirac suggested a 
time-dependent extension (tdHF) that would account for 
dynamical processes Q] which has been the basis of many 
studies of oscillations in systems such as atoms, nuclei 
and molecules. Time-dependent HF theory is not re- 
stricted to small amplitude oscillations, although in this 
limit one obtains a scheme based on ground state pertur- 
bations which is not explicitly time dependent. In fact, 
one can obtain linearized equations (tdHF*) which arc 
essentially equivalent to the well-known random-phase 
approximation (RPA) 

In this Letter we numerically study the tdHF evolu- 
tion for electrons confined in a two-dimensional semicon- 
ductor quantum dot by direct integration of the time- 
dependent equations. We show the equivalence with 
the linearized tdHF* approach previously used by one 
of us and, at the same time, prove the minor role of 
the ground state symmetry-breaking on the far-infrared 
(FIR) absorption. We also explore the nonlinear regime 
by going to the large amplitude limit where we find a 
fast energy transfer from dipole to internal breathing 
and quadrupole oscillations, made possible by the non- 
parabolicity of the confining potential. The excitation of 
internal breathing and quadrupole modes results in size- 
able modifications of the dipole absorption, as compared 
to the small amplitude case. 

As described in many textbooks, the essential ingre- 
dient of tdHF theory is the description of the state of 
the system at any time of the dynamic evolution by a 



single Slater determinant, built with single-particle or- 
bitals. Given initial conditions, the central equations of 
the theory provide the time evolution of the orbital set 
{(pijj{r,t)}. The orbital labels indicate state number (i) 
and spin z-projection {rj =t, J.). The tdHF equations for 
electrons confined to the xy plane by an external v^xti^, t) 
and vector A(r) — B/2{—y,x) potentials read Q 



^(-iV + 7A(ri))Vi;H(ri,i) 



+ t^ea:t(ri,i) + -g*m*^Bsz 



Pr,{Yi,Y2,t) 
dV2— ipiy^{V2,t) , (1) 

where we have defined 7 = e/c and the Hartree potential 
D//(ri,i) = J dr2 p(r2, t)/ri2, given in terms of the total 
density p(r,t). In Eq. (^ we have also defined the one- 
body density-matrix entering the exchange integral 



Pr,(ri,r2,i) = ^ c/?j,,(ri,t)(/3*^(r2,i) 



(2) 



In this work we shall present direct numerical solutions 
of Eq. (|l|) obtained by discretizing the equations both in 
space, into a uniform grid of points, and in time tk = fcAt 
by using the Cranck-Nicholson algorithm [p[ . Specifically, 

(k\ ik) 

by calling and V]^ the local term (within square 
brackets) and the exchange integral in Eq. (Q) at time- 
step k, the algorithm reads 



(fe+i) 



(fc+i) 



+ 



1 

iAt 



lAt, 

h 

2 

(fe) 



(fe) 



2 V 



V. 



(fe) 



(k+l) 



(3) 



(fe+i) 



and V, 



(fe+i) 



As a first guess, we assume /iq 

^i'n which allows an initial solution of (0) . Subsequently 
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are updated and the process repeated 

Sk+i) 



untill convergence in V'i^^^'' is found (typically 4-8 itera- 
tions are enough for a sufficiently small At). 

The system strength function S{uj), or equivalently the 
absorption cross section cr(aj) = ljS{lo) are easily de- 
termined from the time signals (^i'Yi){t) by means of 
a frequency analysis, in a much similar way as used in 
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real-time Density-Functional calculations of atomic clus- 
ters HJ^ and quantum dots jsj. Important differences 
with the latter are the exact treatment of the exchange 
(neglecting correlation though) and the possibility to ex- 
plore large amplitude oscillations since the formalism is 
not centered on a ground-state energy functional. 

In the limit of small oscillation amplitudes we com- 
pare with a frequently used linearized version (tdHF*), 
which provides the linear response to an external oscil- 
lating dipole potential adiabatically switched on 



'(r,t) ^f'^'^Vexp [±ie - i{uj + iT])t] 



(4) 



with (r, 9) the polar coordinates, 77 0+ and '±' la- 
belling the two circular polarization states. The power 
absorbed by the system is 



S{r^(c.)(a|0-|/3)}, 



(5) 



where (j)^'^ is the self-consistent electrostatic potential 
found from the external perturbation (/)°'^* via the dielec- 
tric tensor 101, and 



fO 
J 9 



J a 



(6) 



with /'^ the equilibrium Fermi distribution. Notice that 
within the tdHF* we require the angular part of the or- 
bitals to be simple phases e**^*, which amounts to as- 
suming circular symmetry of the density and mean field. 

Most of the system physical properties are determined 
by the external potential Vextir), which we shall assume 
to be static and circularly symmetric. Following Ref. 

we consider a step in the radial direction separating 
two regions with a different slope in the potential. This 
type of behaviour reproduces to some extent the effective 
potential in quantum dot arrays, where lattice periodicity 
introduces a deviation from a pure quadratic radial law. 
The specific parametrization we use is 

Vext{r) = vi{r) fi{r) + [v2{r) 

+ vi{Ro) - V2{Ro) + A] f2{r) , (7) 

where Rq and A are the step position and height, re- 
spectively. In Eq. (0) vi^2{'r) indicate the potential in the 
r < Rq {r > Rq) region and fi.2(j') is a smooth switching 
function. Namely, these region-dependent functions are 

1 



/i,2M = i flTtanh^— ^ 



(8) 



In the calculations shown below we consider a system 
with N = 6 electrons and Rq = 30 nm, cr = 5 nm. 



uji,2 = 3.37 meV, A = 6 meV. This particular set of 
values is motivated by the experiments in Ref. |9|. 

Figure 1 displays the dipole absorption in the small 
amplitude regime obtained within tdHF and tdHF*. As 
is well known from Kohn's theorem, in purely parabolic 
dots the FIR absorption consists of two bands at energies 



UJ± 



(9) 



where u>c = eB /c is the cyclotron frequency. Looking at 
the lower panel of Fig 1, one can indeed identify a low 
energy peak around 2 meV and a higher-energy structure 
with two peaks. The existence of a structured high en- 
ergy band is in agreement with the experiments and as 
seen in Fig. 1 the two methods yield very similar results, 
showing the equivalence of the very different techniques 
in this limit. Quite remarkably, at _B = IT HF provides a 
symmetry broken ground state, with the six electrons lo- 
calized in a ring-like structure similar to those discussed 
by Yannouleas and Landman ||l0|| or Manninen et al. [0 , 
which is obviously not present within the circular HF*. 
Nevertheless, Fig. 1 shows that the absorption obtained 
from the symmetry-broken ground state does not impor- 
tantly differ from that of the circular solution, proving 
that the details of the FIR absorption are not much sen- 
sitive to the internal structure. For a purely parabolic 
dot there is a complete independence of the FIR absorp- 
tion and the internal symmetry but it is a remarkable 
fact that for the more realistic dot potential (0) the de- 
tails of the high energy branch are not sensitive to the 
internal structure. 

We turn next to large amplitude oscillations within 
tdHF. Figure 2 shows the center of mass trajectory fol- 
lowing an initial rigid displacement to the point indicated 
in the upper right corner. This is a dramatic perturba- 
tion which amounts to an excitation energy of A_B = 26.5 
meV, to be compared with a ground state energy of 
E — 112.7 meV. A relevant observation in Fig. 2 is that 
the amplitude of the oscillation dampens in time, hint- 
ing at an important energy transfer from the center of 
mass to the internal degrees of freedom which are evolv- 
ing in time. To better monitor this coupling mechanism 
we introduce the internal coordinates f = r — Rem, where 
Rem is the center of mass position, and define the inter- 
nal quadrupole and monopole operators as 



XiUk 



ik 

J2 , _2 



N ^ 

ik 



+ ^i^y^ ^^^^ + y^v^ ■ ^^^^ 



Notice that these internal probes involve two-body oper- 
ators when expressed in the fixed laboratory frame. 

Figure 3 shows the time signals corresponding to the 
large amplitude displacement of Fig. 2 as well as the re- 
sults following a much smaller perturbation correspond- 
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ing to the linear regime. As expected from the previ- 
ous figure, the large oscillation dipole signal is reduc- 
ing its amplitude in time, while the internal quadrupole 
and monopole modes get highly excited. The breathing- 
mode result is actually indicating that the dot expands 
after a rather short time, and keeps oscillating around 
the expanded configuration. This energy transfer to 
the breathing mode is not seen in the small amplitude 
regime, where only minor oscillations around the ground 
state configuration are always found. As already men- 
tioned, the mechanism observed here is possible due to 
the non-parabolic confining potential permitting the en- 
ergy transfer between modes, showing that nonlinear pro- 
cesses associated with high energy excitations can effec- 
tively induce strong internal excitations. 

The frequency analysis of the large amplitude signal is 
complicated by the decrease in amplitude. An approxi- 
mate treatment is still possible by defining a time win- 
dow in which the variations are not too large. In fact the 
frequency spectrum evolves smoothly in time, depending 
on the window position. For large times the result is not 
much sensitive on this position and is shown in Fig. 4 for 
B = IT. Notice that in the frequency analysis of the time 
signal each peak is represented by a Lorentzian having an 
artificially chosen width F, for an easier representation of 
the spectrum. Although experimentally rather large F's 
are usually found, in Fig. 4 we show the result obtained 
using either a large or a small value of F, in order to 
clearly display the fine structure of the absorption. 

An important fragmentation appears in the large am- 
plitude spectrum of Fig. 4, due the splitting in many 
closely lying peaks. Actually, the three upper sharp tran- 
sitions of the small amplitude spectrum transform into a 
bunch of peaks. Within a one-particle picture, we at- 
tribute this effect to the large fluctuation of the mean 
field, seen from the internal signals, which causes also im- 
portant variations in the effective^ single particle energies 
and in the corresponding particle-hole transitions. Inter- 
estingly enough, when using a large (and more realistic) F 
we obtain a narrowing of the high energy branch, roughly 
covering the interval [4.75, 5.5] eV. In the small amplitude 
case the higher branch, composed of three peaks, approx- 
imately lies at [4.25, 5.5] eV. The observed effect is thus a 
dynamical narrowing of the high energy branch induced 
by the large amplitude oscillation. 

To summarize, the tdHF equations have been inte- 
grated in time to describe electronic oscillations in two- 
dimensional quantum dots confined by a circular non- 
parabolic potential. In the limit of small amplitudes the 
results are very similar to those obtained in the linearized 
version of the circularly constrained theory (tdHF*). The 
symmetry broken ground states do not significantly mod- 
ify the dipole absorption in this regime. The nonlinear 
behaviour has been investigated by using large amplitude 
initial displacements. An effective transfer from the cen- 
ter of mass dipole to internal monopole and quadrupole 



modes has been observed in this case, thus providing 
a mechanism to strongly excite internal modes in non- 
parabolic dots. Additionally, the large amplitude mo- 
tion has been found to induce a dynamical narrowing of 
the dipole absorption. This effect would manifest in ab- 
sorption experiments using a high-intensity far-infrared 
radiation. 
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FIG. 1. The power absorption for 6 electrons obtained 
within tdHF and using the linear version with circular sym- 
metric tdHF*, for two different values of the magnetic field 
B. 




200 400 600 800 



time (a.u.) 



3 „ 
03 



-1 



-2 



-2 




-1 



X (a.u.) 



FIG. 2. Center of mass trajectory following a largo initial 
displacement to the point in the upper rigth corner, obtained 
within tdHF for the quantum dot considered in this work at 
B = IT. Dotted, thin solid and thick solid lines indicate the 
motion in subsequent intervals of 12 ps, respectively. See 
text. 



FIG. 3. Time dependence of the center of mass (CM) 
dipole, internal quadrupole and internal monopole signals in 
the small (solid) and large (dashed) amplitude regimes at 
B = 1 T. 
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FIG. 4. Dipole absorption corresponding to the small (up- 
per panel) and large (lower panel) amplitude signals of Fig. 

3. Two different peak widths have been used for a bettor pre- 
sentation. The large amplitude results are obtained from the 
time window [500,800] a.u. 
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